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Abstract. Using both theoretical modeling and computer simulations we study a model system for DNA 
interactions in the vicinity of charged membranes. We focus on the polarization of the mobile charges in the 
membranes due to the nearby charged rods (DNA) and the resulting screening of their fields and inter-rod 
interactions. We find, both within a Debye-Hiickel model and in Brownian dynamics simulations, that the 
confinement of the mobile charges to the surface leads to a qualitative reduction in their ability to screen 
the charged rods to the degree that the fields and resulting interactions are not finite-ranged as in systems 
including a bulk salt concentration, but rather decay algebraically and the screening effect is more like an 
effective increase in the multipole moment of the charged rod. 
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^ 1 I ntroduction In this paper we treat this problem using a two-dimensional 

r^ \ salt solution model R,P, pi, follO] to account for the charged 

5t Recent experiments on various systems containing DNA 



strands and charged surfaces have raised an interest in 
understanding how these surfaces screen the electrostatic 



surface (membrane), while the DNA strands are modeled 
as negatively charged rigid rods. We treat charge neu- 
tral systems where the over all charging of the surface 

fields of the DNA and the resulting interactions between , ^, . . ^i r • ^^^ ^^ 

and the rod is zero, thus locusmg either on overall neu- 

strands. The interest emanates from sources varying from ^ „ , , . . i- • . n • . . {TTh ^ • ^ ^ 

trally charged mixed lipid iluid membranes ||12||, or highly 

understanding prokaryotic DNA replication Pl, to non- viral , i r ^ i • i ^i ^ • ^ i 

^^ charged surfaces to which the counter- ions are strongly 

^'t-FQ'tJ' ^ ^yH3' bound and therefore treated within a two-dimensional ge- 
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ometry. We neglect the possible dependence on dielectric 2 Model 

properties of the components in the experimental systems 

The free energy of a system of fixed and mobile charges 
mentioned above in order to simplify the theoretical pic- 
includes electrostatic terms and entropic terms: 
ture. The effects of the thickness and the dielectric prop- 
erties of the layers will be published elsewhere |l3|. p^^^^ = if! / / c^(r)a(rO ^ ^{^)^{^1 ^ ^ (T{r)U{T') \ ^^^^, ^ 

2 e J \ \y-y'\ \y-y'\ \y-y'\ ) 

We focus on very simple geometries in order to under- (1) 
stand how the DNA strand and surface charges interact S = — (cr+(r)(log(cr+(r)ao) — 1) + cr_(r)(log(cr_(r)ao) — l))dY. 
and screen. The geometry is that of one, infinite, salty sur- (2) 
face decorated with either one or two DNA strands. We Here a-^ and cr_ are the number densities of the pos- 
calculate the charge distribution around the DNA and the itive and negative mobile charges where the total mo- 
resulting interaction between the two strands. We calcu- bile charge density is given by ea = ecr+ — ecr_, and 
late the interaction assuming the DNA strands are slightly eU = eU-^ — eU- are the fixed charge densities, 
raised above the surface (See Fig. 1). Minimizing the Grand Potential: 

The theory we use to obtain our analytical results is 



G 



= Felec - ksTS - / /i ((7+ + a_) dY, (3) 



within the Debye-Hiickel approximation] 14 : We minimize 

the free energy of the system with respect to the charge ^^^^ rested to the mobile charge densities, (7+ and a_, 

densities and use the result in the Poisson equation (thus yi^^^Ls. 

obtaining the Poisson-Boltzmann (PB) equation) which , , -, -e0(r)+^ e,p{r)+^l 

^ ^ ^ cr+(r)=ao"^e ^bt cr_(r) = ao"^e^^^. (4) 

we linearize with respect to the electrostatic potential. 

Solving this equation leads to the optimized self consis- ^ere the chemical potential of the positive and negative 
tent charge distribution which we can insert back into the charges, /i, is taken to be equal since we treat the ther- 
free energy in order to obtain the resulting interactions. In modynamic limit of an infinite system where the average 
the next section we introduce the model and the formal number of positive and negative charges are equal far from 
results, and in the following section we apply it to find rods, and 

the interactions between two strands. We compare these e [ f (Jj^(y') — (J-(y') + Ej^(y') — E_(y') 

results with new simulations of two-dimensional salt solu- 
tions and we conclude with a discussion on the limits of 
applicability and relevance of this model. 



w-;/( '-'-''-'-'-;:!gf'-''-"-'-'> -. 

(5) 
is the resulting electrostatic potential. Inserting these re- 
sults into the free energy we find the formal expression: 



G=\j{S{v)-a{v))<t>{v)dv. (6) 
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In this expression we have dropped constant terms that do The potential that solves Eq. ^ has two contributions: a 

Q 



not depend on the exact geometry of the system since we singular part, (I) sing = — 7=^=^5 which solves the equation 



eV^i^ir) = j-{c|)s^ng^^^)S{z). (9) 



are interested in how the energy depends on the distances for a single point charge: eV^0sin^(r) = —4:e7rQ6{p)6{z). 

between the charges objects. We will use this expression The second contribution,?/^, arises from the mobile charges 

in the next section to calculate the interactions in this sys- on the surface and must satisfy the remaining equation: 
tem. Note that the mobile charge density, a, enters with -. 

an opposite sign to what one would have naively guessed 

to be the interaction. This is due to the fact that this Eq. | is the Laplace equation: V^ip = 0, with the special 

term enters as an entropic contribution and therefore in- boundary condition at z = 0: 
dicates how the entropy has been reduced (and thus the 

-.N -. 1 r 1 ^ ^ = \{V^9sing)\z=0' (10) 

tree energy mcreased) due to the arrangement 01 the mo- oz oz A 

bile charges around the fixed charges. ^nr 1 r / -^i r -i r 1 x- r^i t 1 

^ ^ We solve tor ip with a lamily 01 solutions 01 the Laplace 

By inserting the distributions of Eq. || in the Poisson equation- 
equation we get the Poisson Boltzmann equation Jlql: ,, f _n\y\ , X /N 

"^—^ '0(r) = / aqC "^^^Uoiqpjdq. (11) 

/7 -e0(r) + M e0(r) + A^ \ _ \ 

eV (/)(r) = -4e^ |e 'bt _e ^bt J^q -^U{r)jS{z). Using the identity (/),,n^ = --4_ = f / g-^l^ I Jo (^p)dg, 

V' ) the boundary condition (Eq. |l^) is easily satisfied with 

The S function was introduced because the charges are ^^ = -sSli' ^^^ ^^^ ^^^^^ electrostatic potential of the 

confined to the surface at 2: = 0. Linearizing this equation system is given by- 
yields the Debye-Hiickel (DH) equation p5|: 

(l>point = (losing + V^ = Q / e'"^^^^ Jo{qp)^-—-—dq. (12) 

/I \ I 8e^7re^B^ This potential can now be integrated over a line to give 

eV>(r) = -0(r) - 47reZ'(r) S{z) ; 



the potential of a charged rod on a salty surface: 
(8) 



Solving Eq. I for a given fixed charge distribution Z'(r), (j)^^^ = / (j)j,ointMdy = 2r / e"^'^' ^^^''' ' ^ — dq, 

J J q 2e(7A + 1 



_^l^| COs((7x) 2qX 
q 2eqX + 
yields the mobile charge distribution and electrostatic po- (13) 

tentials and fields. In our case the fixed charge is that of where r is the charge per unit length on the bare rod (in 

a uniformly charged stiff rod (model DNA). We first solve the case of DNA: r c^ -e/1.7A). Within this model the 

Eq. I for a fixed point charge Q, {i.e., E{t) = QS{p)S{z)) resulting charge distribution on the surface is found to be: 

and then, since the problem is linear, we integrate to find 



the corresponding solution for an infinitely long rod. ^ttA tt J 2eqX + 1 



"(") = -2^^™'' = -j/ 



cos{qx)dq 
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3 Interactions ulation results which will be described in the following 

chapter, we have to modify the problem slightly so that 
In this section we calculate the effective interactions be- 

the rods are at least slightly raised above the surface and 
tween the components in the system. The surface charges, 

thus do not interfere with the charge distribution around 
which distribute themselves around the fixed charges, screen 

each other. (This is a requirement of the simulated sys- 
to some extent their fields and thus the direct interactions. 

However, the screening is not as effective as that of a three- 
dimensional salt solution where the exponential screening 
leads to a finite ranged interaction. In the case of a two- 
dimensional salt, although reduced, the fields are still long 
ranged H. This can be seen when we take the limit of large 
distances from the rod and calculate the fields resulting 
fromEq. |l|[|0|: 



tem in order to allow for ions to move from one side of the 
rod to the other.) For a rod raised by a small (compared 
with A) distance d above the surface, the amplitudes of the 
modes in i/j are now modified to be aq{d) = —jr^^^x+Vi' 
For rods that are close to each other {D < A) we know that 
the fields will lead to a logarithmic inter-rod interaction. 
However, the fields farther away are more complicated and 

the resulting interaction is not obvious. 

T^ r X ^^ SeVA^ 2erA 

-tjx[x ^ X^ z = O) 0:^ 2 — Ez[x = 0, z ':^ X) c^ — —. The interaction can be calculated numerically as a 

(15) function of distance between the two rods, however, when 

These effective fields are dipolar in nature rather then the the rods are separated by a distance I) > A the po- 

usual 1/r term for a charged line. However, they are not tential (Eq. |l|) and charge distribution, a (Eq. |l^, can 

exponentially screened. (Close to the rod (distances < A) be analytically approximated yielding a simple form for 

the electrostatic potential is not screened and therefore is the interaction energy as a function of the distance D. 

logarithmic as is the case for a bare charged rod.) Once this approximation is made it can be shown that 

(j){x^ z = 0) ^ Xjx^ and g{x) ~ 1/x^ and the integrals in 

3.1 Interaction between two neighboring rods Eq. ^ are simplified to yield the interaction per unit length 

of the rods as a function of distance: 
In order to calculate the interaction energy between two 

G(D X)-^^^^\l -- - T-V I fl6) 

charged rods adsorbed on a salty surface we use Eq. || for ^ ^ ~ ^]j'2 I /^\ 4 yxj ) 

the free energy of the system with the charge distribution Here D is the inter-rod distance and d is the distance 

of two rods separated by a distance D. We make use of between the rods and the charged surface (Fig. 1). 

the fact that the DH equation (Eq. ||) is linear so that we The interaction is similar to that of two rods of dipoles 

can solve for each rod separately and then superimpose at a distance D apart. However, we can not say that the 

the potentials and charge distributions of the combined charge distribution is actually dipolar since the field in the 

system of both rods. In order to compare with the sim- perpendicular direction {Ez) behaves differently. 
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In general one usually expects the interaction between [^ , 

charged objects in solution to be dominated by the os- t\T2 ( f d \ f D 

Urr[x,z) = m ( 2 cosh |27r— 1 — 2 cos ( 27r — 

motic pressure of the solute, in this case the two-dimensional , , 

(18) 



salt solution. However, in this two-dimensional case the 
imperfectly screened electrostatic interactions dominate 
over the osmotic pressure which decays more quickly as a 
function of D ( nosmoUc ^ kBT{cT+{D/2) + a_(I)/2)) c^ 
0{(jP') ^ 0{D~^). Differentiating equation|l6| leads to the 
correct two-dimensional pressure between the rods which 
is dominated by electrostatic contributions: 



and the normalized interaction, Ucr{x, z), between a point 
charge and a rod is given by replacing r2 with q/Ly in the 
above expression, q being the fractional charge of the point 
charge. Energy is here normalized to Eq = e^/47reo. The 
corresponding interaction energy between point particles 
in partially periodic media can be found in Ref. jl^. We 
further employ a short range repulsive interaction poten- 



^r 16f Al^ / S <^ 1 /" r]\^\ ^^^^ '^^^ units of kcal/mol) between ions in the plane: 

UlAt) = < 



^ = -^ = ^D^^nAn a' ' ^''^ 



^"^ L^-^ J (19) 

, otherwise 



4 Numerical Simulations 



with p = 4 and e = 0.01. 

Figure 2 shows the simulation results for two rods of 

charge density, r = e/2A, at a distance d = lOA above 
In order to verify the predictions of the above theory, the surface, and with a distance D between them. The 
we have performed numerical Brownian dynamics Simula- simulations have been performed by initiating the posi- 
tions of the effective interactions between charged stiff (in- tions of the in-plane ions randomly and allowing them to 
finitely long) rods (along the y direction) above (in the z equilibrate for > 10^ time steps {dt = 0.005 in normalized 
direction) a surface, {x^y^z = 0), in which charged mono- time units) before averaging the mean forces over > 10^ 
valent particles can move (see Fig. 1). The simulations steps. We show the attractive mean forces between the 
are performed at a temperature, T = 300K, and with rods and the surface (Fig. 2a) as well as the repulsive 
uniform dielectric constant, e = 80, simulating the con- mean force between the rods (Fig. 2b) for several differ- 
tinuum properties of bulk water. Simulating a non-zero ent ionic strengths of the surface. The prevailing trend is 
salt concentration in the surface, we have applied periodic that the mean force between the rods decays as I)~^ for 
boundary conditions in the plane (x,?/), with a periodic- large D and that the mean force between the rods and 
ity of {Lx^Ly) = (400,40), where length is normalized to the surface approaches the asymptotic (single rod limit) 
To — lA. The normalized long range interactions between as D~^ for large D. This is in agreement with the pre- 
rods in this partially periodic system are thus given by dictions given by Eq. |l^ and ^ It should be noted that 
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given the periodic boundary conditions necessary for sim- charge on the surface); and (A) o and x combined. Even 

ulating a non-zero concentration of salt in the surface, the case where only counter- ions are present shows rea- 

we do not have complete freedom to consider the limit sonable agreement between simulations and theory. We 

D ^ oo. The largest possible D is given by Lx/2 (where have also verified that good comparisons between simula- 

all forces between the rods, in the x direction, are zero by tions and theory hold for other values of d, the distance 

symmetry) and one should therefore only consider simu- between rods and surface. The general trend is that as 

lation results for D somewhat less than Lx/2 in order to d increases, so does the distance, D^ at which the mean 

obtain meaningful results to be compared to the theory, forces approach the predicted slopes shown in figure 2. 
which does not consider a periodic array of particles and 

rods. The results show some signs of insufficient averaging 5 Conclusions 
for large D. Clearly, this is due to the very small mean 



forces that we try to evaluate, combined with the rather 
small number of simulated particles from which the av- 
erages are generated. The small systems are necessitated 
by the long range interactions which require all charged 
objects to interact with all other charged objects, thereby 
increasing the simulation time by the square of the num- 
ber of charged objects simulated. Since the averaging only 
gets linearly better with the number of particles, we are 
limited to small systems if we want to perform calcula- 
tions of the electrostatics carefully. Even so, the simula- 
tions do overall agree with the predictions. In fact, the two 
main predictions, namely that the inter-rod force and the 
rod-surface force asymptotically behave as D~^ and D~^^ 
respectively, seem to be very robust results even though 
they were derived in the thermodynamic limit of a large 
number of overall salt ions in the plane. The three different 
cases shown in figure 2 represent (o) only counter- ions to 
the rods; (x) only salt at concentration 1/400A~^ (charge 
neutrality is satisfied by evenly smearing the rod counter- 



We have presented both theoretical and simulation results 
for the interactions between charged DNA-like rods near a 
salty surface. We have focused on the polarization of the 
surface charge distribution and how it affects the fields 
and resulting interaction between two such rods. Our main 
conclusion is that when the mobile charges are confined, as 
is the case in our treatment, to a two-dimensional surface, 
they do not exponentially screen the fields, and hence the 
effective interaction between the rods is not finite-ranged. 
Both theory and simulations show consistent power law 
interactions both for the force between the rods and for 
the force they apply to the charged surface. 

Despite the fact that the limits of applicability of the 
theory and simulation are not the same: the theory is valid 
at the limit of a large number of ions, and assumes an in- 
finite surface with just two rods, while the simulation is 
restricted, for the reasons mentioned in section Q, to a 
small number of particles and periodic boundary condi- 
tions, we still find a finite region of inter-rod distances 
where the two agree fairly well. This agreement indicates 
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that although the theory is approximate it is robust non Yig. 1. Schematic of the model system of membrane and DNA 

the less. strands 

The systems we treat in this paper may seem relatively ^'S" ^- ^^^'^ ^''"'''^ ^^^^^^^^ '""^^ """"^ ^'^'"'^^'^^ (^) """"^ ^^^^^^^^ 

,.p.,T .11 P1..1 r the two rods (b) as a function of the inter-rod distance, D. 
artificial because the charges are confined to the surface 

. . . ... . 1 1 1 n 1 • 1 System parameters are: Lx = 400, Ly = 40, r = e/1.7A, d — 

and there is no, or very little, residual bulk salt m the sur- 

10, T — SOOiT. Distances are in units of A. Forces in figure (a) 
rounding water solution. In addition, we do not take into 

are relative to the force at D — Lx/2. Results are shown for 
account the effects of the non-zero thickness of lipid mem- 
three different two dimensional salt concentrations (see figure) . 
branes and the relatively low (compared to the surround- 
ing water) dielectric constant of the lipid. However, de- for distances L^ > A, and we do not expect the power law 
spite these limitations, we can apply these results to some forces to change in this regime. 

experimentally investigated systems. Specifically, in re- We are grateful to Cyrus Safinya, Ilya Koltover and 

cent x-ray experiments[| that studied DNA-Cationic lipid Gerard Wang for discussing with us their experimental 

complexes the structures that were found were formed results. We thank Helmut Schiessel, Dov Levine and Nily 

by layers of membranes intercalated by ordered DNA do- Dan for useful discussions. PP and KM acknowledge the 

mains. Most of these experiments were performed in very partial support of the MRL Program of the National Sci- 

low bulk salt concentration, and moreover, because the ence Foundation under Award No. DMR99-72246 in addi- 

counter-ions were trapped between the membranes they tion to Awards 8-442490-21587 and 8-442490-21825. NGJ 

were effectively restricted to two-dimensional space. (The acknowledges support by the Director, Office of Advanced 

lower dimensionality of the space is "measured" relative Scientific Computing Research, Division of Mathemati- 

to the fixed charged objects that polarize the surround- cal. Information, and Computational Sciences of the U.S. 

ing solute. In this case we are studying the screening of Department of Energy under contract number DE-AC03- 

the DNA and therefore the space available to the ions 76SF00098. 
for redistributing themselves is effectively two-dimensional 
compared with the DNA, despite the fact that the ions 
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